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EXTENDING THE RANGE OF LITHOGRAPHIC SIMULATION INTEGRALS 

Cross-Reference To Related Applications 

The present application is related to commonly assigned U.S. patent application Serial No. 

(Attorney Docket No. FIS920030108, titled "INCORPORATION OF A 

PHASE MAP INTO FAST MODEL-BASED OPTICAL PROXIMITY CORRECTION 
SIMULATION KERNELS TO ACCOUNT FOR NEAR AND MID-RANGE FLARE"), 

U.S. patent application Serial No. (Attorney Docket No. FIS920030109, 

titled "IMPROVEMENT OF PERFORMANCE IN SECTOR-BASED OPC ENGINE 
UTILIZING EFFICIENT POLYGON PINNING METHOD AND SIMULTANEOUS 
COMPUTATION OF MULTIPLE SAMPLE POINTS"), U.S. patent application Serial 

No. (Attorney Docket No. FIS9200301 10, titled "RENESTING 

INTERACTION MAP INTO DESIGN FOR EFFICIENT LONG RANGE 

CALCULATIONS"), and U.S. patent application Serial No. (Attorney 

Docket No. FIS920030262, titled "SIMULTANEOUS COMPUTATION OF MULTIPLE 
POINTS ON ONE OR MULTIPLE CUT LINES"), filed on even date herewith, which are 
hereby incorporated by reference herein in their entirety. 

Background Of The Invention 

1. Field of the Invention 

This invention relates generally to the field of optical microlithography, and particularly 
to the use of triangle convolution to a sector-based OPC engine. More particularly, the 
invention relates to the use of triangle convolution on unbounded sectors, which 
inherently regularize the convolution of an intermediate-range, long-range, or infinitely 
extending kernel without spatially truncating the kernel, nor requiring excessive ROI 
size. 

2. Description of Related Art 

The optical microlithography process in semiconductor fabrication, also known as the 
photolithography process, consists of duplicating desired circuit patterns as best as 
possible onto a semiconductor wafer. The desired circuit patterns are typically represented 
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as opaque and transparent regions on a template commonly called a photomask. In optical 
microlithography, patterns on the photomask template are projected onto photoresist- 
coated wafers by way of optical imaging through an exposure system. 

Aerial image simulators, which compute the images generated by optical projection 
systems, have proven to be a valuable tool to analyze and improve the state-of-the-art in 
optical lithography for integrated circuit fabrication. These simulations have found 
application in advanced mask design, such as phase shifting mask (PSM) design, optical 
proximity correction (OPC), and in the design of projection optics. Modeling aerial images 
is a crucial component of semiconductor manufacturing. Since present lithographic tools 
employ partially coherent illumination, such modehng is computationally intensive for all 
but elementary patterns. The aerial image produced by the mask, i.e., the light intensity in 
an optical projection system's image plane, is a critically important quantity in 
microlithography for governing how well a developed photoresist structure replicates a 
mask design. 

Mathematical difficulties arise in optical proximity correction calculations that, in part, 
can be attributed to divergence issues when attempting to analytically express imaging 
kernels over the unbound sectors into which mask polygons are conventionally 
decomposed. Attempts to account for incoherent flare are also subject to these analytical 
restrictions. 

Accuracy is of critical importance in the computation of calibrated optical or resist 
models. The accuracy in the simulation of wafer shapes is necessary for a better 
understanding and evaluation of the correction methodologies. Through analytical 
processes, fidelity in the wafer shapes to the "as intended" shapes may ultimately achieve 
better correction of the mask shapes. An increase in yield during chip manufacturing is a 
direct consequence of achieving this accuracy. 

The present invention introduces the utilization of triangle convolution to a sector-based 
optical proximity correction engine. Using triangles in the mathematical model to 
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eliminate the difficulties of unbounded comers allows for integrals that acconmiodate 
intermediate-range and long-range effects. 

In the prior art, the following mathematical treatment in the optical proximity correction 
engine is commonly used. These approaches are all, in one form or another, related to the 
Sum of Coherent Source (SOCS) method, which is an algorithm for efficient calculation 
ofthe bilinear transform. 

Sum of Coherent Source (SOCS) method 

The image intensity is given by the partially coherent Hopkin's equation (a bilinear 
transform): 

lo (r) = JjJJdr'dr''h(f -rOh-(f-Oj(r^ rV(r^^^ 

where, 

h is the lens point spread function (PSF); 
j is the coherence; 
m is the mask; and 
lo is the aerial image. 

By using the SOCS technique, an optimal n-th order coherent approximation to the 
partially coherent Hopkin's equation can be expressed as 

Io(7) = Z^.h®A)W| 

Where /l^,^^(r) represents the eigenvalues and eigenvectors derived from the Mercer 
expansion of: 

WCF',^ = h(F)h-(Oj(F-0 = Z^A(r')*k(n 

k=l 

which suggests that a partially coherent imaging problem can be optimally approximated 
by a finite summation of coherent imaging, such as linear convolution. 
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SOCS with pupil phase error 

The above calculation assumes an ideal imaging system. However, when lens aberration 
is present, such as pupil phase error and apodization, one must include the pupil function: 



Because of the possible higher spatial frequency in the wavefront function, h(f) will 
have a larger spatial extent. In this case, the number of eigenvalues and eigenvectors 
required will be higher than the ideal system. Hence, the kernel support area is extended 
to take account of the contribution from a distance greater than XfNA, However, the basic 
mathematical structure and algorithm remains the same. 

Physical model of Flare 

Flare is generally described as the image component generated by high frequency phase 
"ripples" in the wavefront. Flare thus arises when light is forward scattered by appreciable 
angles due to phase irregularities in a lens. Such irregularities are often neglected for three 
reasons. First, the wavefront data is sometimes taken with a low-resolution interferometer, 

and moreover may be reconstructed using an algorithm of even lower resolution. Second, 
even when the power spectrum of the wavefront is known or inferred, it is not possible to 
include the effect of high frequency wavefront components on an image integral that is 
truncated at a short ROI distance, causing most of the scattered light to be neglected. Last, 
it is not straightforward to include these terms in the calculated image. The present 
invention addresses these problems. 

Bearing in mind the problems and deficiencies of the prior art, it is therefore an object of 
the present invention to provide a method for calculating intermediate and long-range 
image contributions from mask polygons. 




where. 



?{p) is the pupil transmission function; and 

W(/5) is the pupil phase function, which contains both aberration and 



flare information. 
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It is another object of the present invention to apply triangle convolution techniques to a 
sector-based OPC engine. 

It is a further object of the present invention to provide a method to calculate incoherent 
flare. 

It is yet another object of the present invention to reduce the time needed for numerical 
integration for intermediate range calculation. 

Another object of the present invention is to provide an analytical solution for the triangle 
convolution on a power law incoherent kernel. 

A further object of the present invention is to present a programmable method to divide 
long-range and intermediate-range calculations by truncating mask but not kernel. 

Another object of the present invention is to provide a method to account for the image 
component generated by high frequency phase "ripples" in the wavefront. 

It is yet another object of the present invention to provide a method to calculate the 
effects of phase irregularities in a lens. 

Another object of the present invention is to provide a method that includes the effect of 
high frequency wavefront components on an image integral that is truncated at a short 
ROI distance. 

Still other objects and advantages of the invention will in part be obvious and will in part 
be apparent from the specification. 

Summary of the Invention 

The above and other objects, which will be apparent to those skilled in art, are achieved 
in the present invention, which is directed to in a fust aspect, a method of performing 
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model-based optical proximity correction comprising: calculating a finite integral for 
each sector of a mask polygon having defined comers to determine sector contributions 
from vertices of said mask polygon; integrating on a finite shape in said sector; simiming 
up the contribution of said sectors of said mask polygon; and calculating interaction 
effects from said finite integral calculation. The finite shape may comprise at least two 
triangles. Each of said comers of said polygon is integrated. The method includes using a 
triangle convolution on unbounded sectors or a triangle convolution on a power-law 
incoherent kernel. The triangle convolution is reduced to an analytic expression that 
contains an incomplete beta fiinction to facilitate numerical integration. A minimum 
radius, Rmin, is determined such that said convolution does not double count any 
contribution from inside a ROI. 

The method further includes providing a mle for a positive or negative sign of said 
contributions, comprising: assigning said sign of a first triangle to +1 for two partitions 
on a side of a first edge that are inside a sector and for said partition opposites, otherwise 
assigning said sign of said first triangle to -1; assigning said sign of a second triangle to 
said sign of said first triangle if said first and second triangle does not overlap, otherwise 
assigning said sign of said second triangle to the opposite of said sign of said first 
triangle. A hole is analytically defined in a flare kemel having said radius Rmin and any 
contributions from said flare that are within Rmin are treated with a partial coherent 
bilinear convolution. The radius Rmin is approximately equal to ^NA. Calculating said 
interaction effects includes integrating over a coarse grid of a plurality of mask polygons. 

In a second aspect, the present invention is directed to method of extending lithographic 
simulation integrals to include intermediate-range and long-range distance scales, said 
method comprising: dividing and decomposing at least two mask polygons into sectors 
made up of triangles having comers and edges; performing a distance scale calculation 
from comers of a plurality of said at least two mask polygons; implementing a polygon 
pinning algorithm; performing a triangle convolution of said plurality of mask polygons; 
combining contributions of said sectors to yield a convolution of said plurality of said 
mask polygons; combining said plurality of said mask polygons; and computing an aerial 
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image for optical proximity correction from said combined sector contributions. The step 
of performing a distance scale calculation includes performing finite integrals from semi- 
infinite comers of said mask polygon. The method ftuther comprises: calculating 
contributions from said each of said polygons at an image sample point; converting area 
convolution into a position vector dependent integral around a feature perimeter using 
first and second position vectors; assigning a resultant indefinite integral for said edges, 
such that a contribution from an edge between vertices of said fu*st and second position 
vectors is a difference of indefinite integrals evaluated at each position vector; and 
parameterizing said resultant indefinite integral as function of an orientation of said edge, 
a perpendicular distance from said position vector to said edge, and a polar angle defined 
from said first position vector to said second position vector. A rule for a positive or 
negative sign of said contributions is also provided. The rule is adapted for said sectors 
including six basic sectors comprising 45-degree and 90-degree sectors. Values of said 
convolution are stored in a first look-up table for various values of said position vectors. 
The rule comprises: assigning said sign of a first triangle to +1 for two partitions on a 
side of a first edge that are inside a sector and for said partition opposites, otherwise 
assigning said sign of said first triangle to -1; assigning said sign of a second triangle to 
said sign of said first triangle if said first and second triangle does not overlap, otherwise 
assigning said sign of said second triangle to the opposite of said sign of said first 
triangle. The method ftirther includes storing said signs of said triangle convolutions 
within a second look-up table, said second table having a row index and a colunm index, 
such that said row index is assigned for said first edge of said sector and said column 
index assigned for a second edge direction, uniquely specifying said sector type, said row 
and column indices determined for any sector in consideration, based on said first and 
second edge directions. First and second sign matrices, Mat 1 (edge l,edge2) and 
Mat2(edgel,edge2), are provided for said second look-up table, said matrices expressed 
as a function of said first and second edges, such that said first sign matrix. 
Mat 1 (edge l,edge2), determines which of said partitions where said direction vector may 
lie will have a positive sign for said first triangle, and said second matrix, 
Mat2(edgel,edge2), determines which of said partitions where said direction vector may 
lie gives the same sign for said second triangle compared to said first triangle. The row 
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and column index of each matrix are given by eight basic partitions comprising: up, 
down, right, left, up-right, up-left, down-right, and down-left. The method fiirther 
comprises multiplying said sign by -1 if said sector is convex, or by +1 if said sector is 
concave. The first matrix signifies which subsets of said partitions carry a positive sign or 
negative sign, such that if said partition matches a matrix entry, said sign of said first 
triangle convolution is assigned a positive value. The sign of the second triangle 
convolution is set equal to said sign of said first triangle partition within which said 
direction vector resided, else said sign is set opposite said first triangle convolution sign. 
A constant of integration is added in an indefinite integral for each edge, such that said 
constant of integration is independent of said polar angle, said constant of integration 
allowing for removal of any infinities in contributions from vertices in isolation. The 
constant of integration is present only in said indefinite integral, vanishing when a 
difference is taken of said indefinite integral at two endpoints of said edge. A 
contribution from an individual mask polygon that has a plurality of edges is represented 
by a simimation over each edge of said mask polygon such that for each comer of said 
edge said constant of integration is subtracted from said indefinite integral. Singularities 
that arise fi-om unbounded individual comers are resolved by adding constant of 
integration terms that leave comer contributions fuiite, such that each comer contribution 
is a summed flare from a pair of said triangles of said triangle convolution. The comer 
contribution comprises a calculation of a contribution from two of said triangles per 
vertex, such that said contributions from said triangle convolution gives a net 
contribution of said mask polygon after summation through all vertices. 

In a third aspect, the present invention is directed to a program storage device readable by 
a machine, tangibly embodying a program of instructions executable by the machine to 
perform method steps for performing model-based optical proximity correction, said 
method steps comprising: calculating a fmite integral for each sector of a mask polygon 
having defined comers to determine sector contributions from vertices of said mask 
polygon; integrating on a finite shape in said sector; summing up the contribution of said 
sectors of said mask polygon; and calculating interaction effects from said finite integral 
calculation. 
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Brief Description of the Drawings 

The features of the invention believed to be novel and the elements characteristic of the 
invention are set forth with particularity in the appended claims. The figures are for 
illustration purposes only and are not drawn to scale. The invention itself, however, both 
as to organization and method of operation, may best be understood by reference to the 
detailed description which follows taken in conjunction with the accompanying drawings 
in which: 

Fig. 1 depicts a schematic of the scattering of Hght due to a high frequency wavefront 
aberration. 

Fig. 2 depicts a flowchart of the proposed algorithm of the present invention for 
extending the range of lithographic simulation integrals to include intermediate and long- 
range distance scales. 

Fig. 3 depicts the geometry for calculating each comer contribution for the summed flare 
value from a pair of triangles. 

Fig. 4 illustrates the geometry for the derivation of the triangle convolution under the 
condition h > Rmin- 

Fig. 5 illustrates the geometry for the derivation of the triangle convolution under the 
condition h < Rmin. 

Fig. 6 illustrates the geometry for the derivation of the triangle convolution xmder the 
condition h/cos 9 < Rmin- 

Fig. 7 is a schematic of the twenty-four different sectors depicted in six basic groups for 
sector convolution. 



FIS920030107US1 



9 



Fig. 8A depicts two overlapping triangles having opposite signs for use in sector 
convolution. 

Fig. 8B depicts two non-overlapping triangles having the same sign for use in sector 
convolution. 

Fig. 9A depicts an example of sign generation for an orthogonal (90 degree) vertex. 

Fig. 9B depicts and example of sign generation for a 45 degree vertex. 

Fig. lOA graphically illustrates the eight partitions of a mask space for feature edges of 
90 and 45 degree orientation. 

Fig. lOB graphically illustrates a sector defined by two edges following the edge 
convention of the present invention. 

Fig. IIA illustrates the long-range convolution performed on a coarse grid, with each 
grid point representing an average pattern density of patterns inside that particular grid 
point. 

I 

Fig. 1 IB depicts the near range convolution performed via polygon convolution inside a 
domain. 

Fig. 12 is a schematic depicting the spatial extent of the short range and mid-range 
calculation. 

Description of the Preferred Embodiment(s) 

In describing the preferred embodiment(s) of the present invention, reference will be 
made herein to Figs. 1-12 of the drawings in which like numerals refer to like features of 
the invention. 
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This disclosiire describes a new method for calculating long-range image contributions 
from mask polygons. The method works even with imaging kernels that are divergent over 
the unbounded sectors into which mask polygons are conventionally decomposed. It can 
therefore be used to calculate incoherent flare of the kind modeled by a power-law; the 
envelope of a fractal power spectrum. Alternatively, the integrated quantity can be a 
deterministic point spread function (PSF) calculated from a measured wavefront. 

The method, however, is not limited to the integration of optical flare intensity. For 
example, image amplitude contributions may be calculated by convolving mask polygons 
with an eigenkemel from a bilinear integral, such as the Hopkin's Integral that governs 
partial coherent imaging. The kernel may also be non-optical in origin, and empirical in 
nature. It may represent a phenomenological linearization of chemical interactions 
between different printed features in a resist. It is demonstrated below that both optical and 
non-optical effects may be handled in the same formaUsm, so that the present invention 
may be effectively applied to either case. 

Modifications to the Sum of Coherent Source (SOCS) method are proposed to incorporate 

a resist kernel, requiring only a sum of linear convolutions, or alternatively, incorporating 
an incoherent kernel. 

SOCS with resist kernel 

The resist image may be defined by convolution of Io(r) with a resist kernel K: 



It can be shown that the above equation can be formulated as a SOCS decomposition, 



I(R) = fJdfIo(f)K(R-r) 

= j| jj Jj dr dr ' dr" K(R - f ) h(r - ?') h* ( r - r") j(r - r") m(r') m* (r ") . 



(1) 



yielding: 




(2) 



where, 



FIS920030107US1 



11 



u^(R) = J|dfm(f-R)i|/j(f). 



(3) 



gj and \|/j(f) are the eigenvalues and eigenvectors obtained from: 

8jN'i(^') = -Ifd^V(?'^')j(^-?')Vj(r). (4) 

where, 

V(f ,f ') = j|dR'K(R') h(f . RO h*(f ' - R^) . (5) 

This means that with the resist kernel, the final image intensity can be treated similarly to 
the SOCS method, which requires only a sum of linear convolutions. 



SOCS Plus Incoherent Kernel 
By defining 

r' + r" ^ 



r = r - r 



Avg^ 



and using the approximation 



hir - •^)h * (r + — ) = h{P) at large values of r = r - r^^ , 



the Hopkin's Integral becomes: 



ROI 



+ 



outside 
ROI 



(6) 



(7) 



outside 
ROI 



(8) 

where, 

is the component coherent kernel. 

This approximation has a simple interpretation. The image contribution from mask 
polygons within the ROI is accounted for using the usual formulation for partial coherent 
imaging (in the SOCS expansion), whereas contributions from polygons outside the ROI 
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are accounted for in a more approximate way. The contribution from every mask point is 
thus accounted for with different degrees of approximation without double counting. 
Though the integration method changes abruptly at the ROI boundary, it is an 
improvement over the present practice in which contributions from outside the ROI are 
dropped to zero. 

The second integral is not truly an incoherent integral, since the mask-like term in square 
brackets depends on the coherence y(A). An incoherent case corresponds to delta-function 
coherence, which ultimately is not a good approximation when the mask consists of low-k 
lines. In lithography, the total source energy is usually normalized to 1, which means that 
with finite coherence and a chrome mask the quantity in square brackets will always be 
less than the pure incoherent case. When the incoherent assumption is sufficiently 
accurate, the second integral is a linear convolution, which is well suited for the method of 
the present invention. 

The present invention is described in three distinct parts. First, an algorithm is introduced 
that applies to Optical Proximity Correction in optical lithography. A finite integral for 
each sector of polygon replaces an infinite integral used in the prior art that otherwise 
diverges for specific user-defined OPC kernels. Integrating over two triangles, rather than 
integrating on the fiiU sector, achieves a finite integral. Even if the kernel is extremely 
long range, the integral still will not diverge because the finite area is integrated. Second, 
an analytical approach is presented for a power law kernel that reduces the numerical 
integration of a sector to an analytical expression evaluation, which greatly enhances 
computation speed. Third, a methodology is introduced to divide the mask into regions to 
calculate interaction effects, such as intermediate-range and long-range effects, by 
truncating the mask instead of truncating the kernel function as suggested in prior art. This 
division efficiently facilitates the calculation of long-range and intermediate-range 
proximity effect calculations and tailors well to an OPC environment application. 

Different scales of wavefront spatial frequencies 

Fig, 1 depicts a schematic of the scattering of light from a high frequency wavefront 
aberration. High frequency wavefront aberrations 10 are depicted for m cycles of 
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wavelength such that the scattering 12 is a function of the wavelength, the number of 
cycles, the frequency, and the numerical aperture. A harmonic component of phase 
modulation in the pupil that has periodicity m cycles across the nxmierical aperture (pupil 
radius) will deviate light by a scattering distance mX/NA in the image plane. Thus, if the 
radius of the optical kernel is set at 0.5)im, the neglected intensity will include mid-range 
spatial frequencies as well as long-range flare. Indeed, the contributions from all phase 
frequencies much beyond the first thirty-seven Zemike terms will be neglected. 

These neglected frequencies are relevant to diagnostic measurements that concern the 
impact of large pads on the widths of nearby fine lines, generally at separation distances of 
approximately Ijiim or more. Furthermore, it may be desirable to improve the accuracy of 
the primary partial coherent image calculation by modestly extending the range of 
integration. Considering slightly larger scales, it is generally advantageous to model the 
kind of printing differences that are observed when fill patterns are placed around test 
structures at distances beyond approximately 2|im. At the other extreme, non-negligible 
flare contributions are known to arise over a range of many millimeters. 

Generally, power-law exponents of small absolute value characterize data across long 
scales, while terms with larger exponents contribute across shorter ranges. Lithography 
lenses provide large Strehl ratios where very low frequencies (aberrations) are concerned. 
This implies that only small amounts of light can be deviated by features which are within 
approximately XfNA of the image point of interest. Thus, one would not expect to see 
significant flare levels that are governed by power law exponents which approach three, 
because appreciable flare intensities can only accumulate over distances that are large 
compared with X/NA. In general, power law exponents on the order of 2.5 govern 
incoherent flare. Short-range wavefront ripple, beyond the standard thirty-seven Zemike 
terms, may still be important since it affects the coherent or partially coherent interaction 
between nearby features, particularly when the source provides a narrowly structured 
coherence, such as a low sigma, aggressive annular quadrupole. It is usefiil to adopt 
empirical kernels that take into account other medium-range interaction effects, which 
may be non-optical in nature, such as resist loading effects. 
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Typically, it is more efficient to integrate contributions over very long ranges using a grid 
that is much coarser than typical mask polygons. This permits an efficiency improvement 
because it is possible to carry out convolutions with a pattern that is locally averaged on a 
periodic grid (density-map) for many image points at once. Fast Fourier Transforms 
(FFTs) are the conventional method for doing fast bitmap convolutions. Direct integration 
with mask polygons is appropriately carried out in an intermediate-range. 

Non-flare Effect 

Important non-flare effects arise from resist processing. The origin of the resist effect may 
be from photoacid/base diffusion, developer loading effect, or airborne contamination. 
The resist effect is usually modeled by a fiill resist chemical model or a lumped parameter 
model. As an approximation, the resist effect may also be described by an incoherent 
convolution of the image intensity with a resist kemel, experimentally calibrated. One can 
perform a linear convolution of the resist kemel on the calculated image intensity, or 
absorb the resist kemel in the bilinear transform and do a sum of coherence source (SOCS) 
decomposition and perform convolution between the mask function and coherent 
eigerikemels as mentioned above. 

Fig. 2 depicts a flowchart of the proposed algorithm 20 of the present invention for 
extending the range of lithographic simulation integrals to include intermediate and long- 
range distance scales. The mask is divided for long-range calculations 22 and intermediate 
range calculations 24. Polygon pinning 26 is performed under the intermediate-range 
calculation, decomposing the polygon sectors 28, and performing triangle convolution 30. 
A description of a polygon piiming algorithm is defined and taught in a related, cross- 
referenced application, identified above, and incorporated herein by reference. Sector 
contributions are then combined 32, and the aerial image is computed for optical 
proximity correction 34. 

Finite Integrals fi'om Semi-Infinite Comers 

A key component of the present invention is to obtain finite integrals from semi-infinite 
comers. Sector contributions fi-om the vertices of mask polygons are regularized, with tiie 
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polygon convolved with a kernel that is divergent over the sector, such as a power-law 
kernel r"^ with y ^ 2. It is known in the lithography art that mask polygons may be 
represented as a superposition of open sectors having essentially semi-infinite comers, but 
that these sectors are unbounded in some directions, potentially leading to infinite 
contributions. 

Of course, the individual sectors are only a means to an end, namely finding the 
convolution over a finite polygon. Thus, one solution is to bound all sectors at a distant 
perimeter, for example, 301 mm for wafers of the current state of the art. However, this 
technique impacts accuracy. 

A preferred solution is based on the following considerations. Given an observation point 
Xq at which to calculate the contribution from a polygon, a Stoke's or Green's theorem is 

employed to covert the area convolution into an Xq - dependent integral around the feature 
perimeter, which is then the sum of integrals along each edge. 

Assigning the resulting indefinite integral for some edge as Q, the contribution fi-om the 
particular edge between vertices x„ and 3c„^, is then The term may be 

parameterized in terms of three quantities, namely: v|/(x^^, -x^), which is the orientation 
of the edge; /i(^05^n>^ii+i)» which is the perpendicular distance fi*om x^ to the edge; and 
Q(^oj^ii'^fl+p^ii+i)' which is the polar angle fi-om Xq to 3c„^j, with the perpendicular 
defining the 0 = 0 origin. These quantities are illustrated below. Note that 0 is the only 
argument of Q that is dependent on the specific location of the individual endpoints along 
the edge. 

One may also include a constant of integration Q in the indefinite integral for each edge, 
provided Q is independent of 0. The constant of integration is only present in the 
indefinite integral, vanishing when the difference is taken of the values of the indefinite 
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integral at the two endpoints of the edge. Even so, this constant of integration allows for 
removal of any infinities in the contribution from vertices in isolation. 

Adding a constant of integration to the contribution from each edge, the integral over the 
polygon is expressed as a sum over edges: 

n edges 

(9) 

Here, F represents the contribution from an individual mask polygon, for example from 
flare, that has n edges, with each of the bracketed terms in Eq. (9) representing the 
contribution from one endpoint of the edge. F can be rearranged into a sum over the 
vertices by regrouping the terms. In this way, a sum over unbounded sectors is achieved, 
as follows: 

n comers 

(10) 

Note that even though the constant of integration no longer cancels within each term in 
this new summation form (sum over comers), it still cancels over the polygon as a whole. 
In fact any integration constant may be employed as long as it obeys the above functional 
dependence. This is key to resolving the singularity that arises from unbounded individual 
comers. Even though Q„^„_i - Q„+i^„ may be infinite, the sum can be transformed by adding 

Q terms that leave each comer contribution finite, without changing the total. 

The easiest way to accomplish this is to make each comer contribution the summed flare 
from a pair of triangles. The basis for this strategy is the fact that there is one value of 6 
that will always be independent of the endpoints, namely 9 = 0, since 0 = 0 is 
unambiguously defined by the perpendicular from the edge to the point of interest. One 
may then choose Q{^f„,„.^,K^^_^) =Q (M/^^„.p/7„„.,,0), and with this choice, the square- 
bracketed terms are no longer unbounded. Fig. 3 depicts the geometry for calculating each 
comer contribution for the summed flare value from a pair of triangles 32, 34. 
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The derivation of the triangle convolution is identified below. Referring to Fig. 4, and 

= G(r,^),if r^lL^ 

assuming the kernel ^" (11) 

= 0, if r<R^ 

let L be defined as the integral for triangle ACDA depicted as the crosshatched area 42, 
with h and 9 as shown. For the case of h > Rmin (case 1) : 

e riff) 

L = l I G(r,9)rdrd9 (12) 

0 Amin 



where cos 6 ~ — ^ (13) 



For the case of h < Rmin, (case 2) as shown in Fig. 5, the triangle integral is given by 
L\K9) = Triangle ACDA Triangle CDEC. 

LXh,0) = j I G{r,e)rdrde 

=L{h,e)-L{h,e,) 

where cos Oq = ( 1 5) 

For the case where — - — < R min, (case 3) as shown in Fig. 6, the triangle integral is 
cos^ 

given by L"{h,6) = 0 , since the whole triangle is inside the hole of the kernel. 



Rules for sector convolution 

The polygon convolution method involves the calculation of the contribution from two 
triangles per vertex as previously mentioned and illustrated in Fig. 3. The triangle 
contributions give the net of the contribution of the polygon after summation through all 
vertices. Since the polygon can be of arbitrary shape and the observation point Xo can be 
inside or outside the polygon, a rule for the sign of the contribution of each triangle is 
required. In most cases, lithographic mask patterns contain polygons formed by 
horizontal, vertical, and 45-degree line segments only. Therefore, an algorithm must 
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handle 45-degree sectors in addition to 90-degree sectors. This makes for 96 possible 
types of sectors that may exist in any polygon. By taking into account the symmetry 
between the clockwise and counterclockwise vertex transversal direction and the anti- 
symmetry between convex and concave sectors, the unique number of sector reduces to 
twenty-four. The sector is convex when the inside angle of the vertex is less than or equal 
to 180 degrees; otherwise, the sector is concave. Moreover, further study shows that the 
sector convolution is identical in magnitude among four other specific sectors. The 
twenty-four different sectors are depicted in six basic groups A-F, as shown in Fig. 7. 

The numerical value of the convolution, which is always non-negative, is identical for the 
sectors within each group since the triangles created are identical, reducing the number of 
unique sectors to the six basic sectors. Once the sector convolution has been evaluated for 
these six basic sectors, it can then be stored in look up tables for various Xq. Note that any 
one of the sectors in each group can become a basic sector as along as a corresponding 
rule for the sign convention is built. The sign, of these triangle convolutions can be 
depicted within a lookup table also. For 45-degree sectors, extra complexity arises for the 
sign convention of a triangle of the pair. Figure 8 shows two different scenarios for the 
sign assignment. The signs of the two triangles can be the same or opposite, depending 
on whether the two triangles overlap as shown in Fig. 8A, or remain adjacent but not 
overlapping, as shown in Fig. 8B, which in turn depends on where Xq is located with 
respect to the sector. By generalizing this extra condition, a set of sign assignment tables 
can be built to determine the sign of the two triangles. 

General Rule for Sign Determination 

Referring to Figs, 9A and 9B, the general rule for sign determination for both 90 degree 
and 45 degree sectors is ascertained as follows. If 4^1 is assigned to the sign of the first 
triangle, and ^2 assigned to the sign of the second triangle, the rule for Tl is simply 
stated as Tl = +1 for the two partitions on the side of Edge #1 that are inside the sector, 
as well their opposites. Otherwise ^1 is assigned to -1. If the sector is concave, ^1 is 
multiplied by -1, otherwise there is no change as it is multiplied by +1. A sector is 
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concave if its inside angle is greater than 180 degrees, or else it is convex. The rule for 
^2 is simply 4^2= 4^1 if the two triangles do not overlap; otherwise 4^2= -4^1. 

The sign determination process is determined in the following manner. First, two sign 
matrices for table lookup are generated. The row and column index of each matrix are 
given by eight basic directions: up, down, right, left, up-right, up-left, down-right, and 
down-left. The row index is assigned for a 1^^-edge direction of the sector and the column 
index assigned for a 2"^-edge direction. The row and column indices uniquely specify the 
sector type. The signs for the triangle convolution are stored in the matrix when Xo is in 
each of the eight partitions in real space. The first sign matrix elements, generally 
referred to as Mat 1 (edge l,edge2), determine which of the eight partitions that Xo may lie 
will have a positive sign for the first triangle. The second matrix, generally referred to as 
Mat2(edgel,edge2), determines which of the eight partitions where Xo may lie gives the 
same sign for the second triangle compared to the first. 

Next, the row and column indices are determined for any sector in consideration, based 
on the direction of the first and second edge of the sector. The first matrix signifies which 
subsets of the eight partitions carry a positive sign or negative sign. The partition that Xq 
resides in is determined for a particular Xo. If this partition number matches the matrix 
entry, then the sign of the first triangle convolution will be positive. 

If the sector is convex (or concave), the resultant sign is multiplied by -1 (or +1) 
according to the convention used. 

The same row and column indices are then used in the second sign matrix, Mat2. If the 
partition within which Xo resided matches the second matrix entry, the sign of the second 
triangle convolution is same as first triangle, otherwise it is the opposite sign. 

Fig. 10 graphically illustrates the procedure to find the sign of an individual triangle 
convolution. Fig. lOA shows the eight partitions 1-8 of the mask space for feature edges 
of 90 and 45 degree orientation. Fig. lOB shows as an example a sector defined by two 
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edges. The row index is in the up-right direction 102, and the column index in the right 
direction 104. In this example, the aforementioned procedural steps require that regions 
2, 3, 6, and 7 have a positive sign, and the remaining regions negative. Since x© is in 
region 1 in this example, it is assigned a negative sign. However, if this sector were 
concave, the sign would be opposite. The sign remains unchanged when the sector is 
convex. According to the second sign table, region 1 corresponds to a positive sign, thus 
the sign of the second triangle will be the same as the first. 

Polvgon Convolution Using Power Law Kernel 

As mentioned above, the flare effect may be treated as an incoherent kernel as long as 
one is considering contributions outside the coherence distance, Rmin> of the point spread 
function. The incoherent kemel for the flare calculation generally obeys a power law of 
the form G(r): 



Evaluation of the convolution at an observation point Xq for a general kemel usually 
requires numerical integration. However, for a power law kemel, the triangle convolution 
Q(^n,n-i, hn,n-i. dn,n-i) cau, after some algebraic manipulation, be reduced to an analytic 
expression that contains the incomplete beta function, denoted Ix (a,b). Beta is a constant 
evaluated from a beta function. The triangle convolution is represented, for case 1, 2, and 
3 mentioned in the previous section above, as 



G{r)^—whenr>R^ 



(16) 




(17) 



QQ¥,h, ^) = ^i:,(A,^),i = 1,2,3 



(18) 



Incase 1, 



FIS920030107US1 



21 



e ri0) 



0/?min ' 

^ d0 



ri0) 



^min 



(19) 



so, 



where cos^ = - 



r-2 



2-r 



2-r _ u^-r 



9 



2-r 

y-\ 1, 



^1 



2-/ 



/2(r-2) 



(20) 



(21) 



In case 2, 



e r(e) 



L^{h,e)= J \\rdrde 

tfg Rmm ^ 



In case 3, 



1 r-1 1 y-\ 



>/2(/-2) 
(22) 



Z,(^,/«,^) = 0 



(23) 



The above expressions are only valid for j*!. For y=2, the function Li reduces to the 
following, 

L,( h. ^) = ^ \e hi(2h/R„ J -1/2 %n-ie)\ (24) 
where f(x) is the Clausen Integral. 
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The incomplete Beta function and the Clausen Integral may be predetermined for a 
specific y value, and stored for table lookup during any subsequent calculation of Q. The 
analytical expression for Q speeds up the convolution calculation significantly with the 
look up table, avoiding time-consuming numerical integration. The incomplete Beta 
function is derived from the following expression: 
e 



Let = cos (p, dt = -2cos^sin^d(p 

cos-0 Ay-2)I1 



j= f ' 

I -2(1-0'"^" 

J = -Beta{^,-) 
2 2 2 



-dt 



■'Wtf*' 2 '2 2 '2 



where Beta is the Beta function 

1 

Beta{z,w) = ^r-\\-tY~'dt and 

0 

Ix is the incomplete Beta function with argument x: 
I,{z,w) = ]r\\-tr'dt 



(25) 



(26) 



(27) 



(28) 



Using the identity: 



and equation (26), yields: 



I,{a,b) = \-I,.,{b,a) 



(29) 



J = -Beta(^,-) 
2 2 2 



y-l 1 



(30) 



Effect of Rm in 

Rmin is determined such that the polygon convolution does not double count the 
contribution from inside the ROI, which essentially corresponds to the diffraction-limited 
PSF and the first 37 Zemike terms of aberration. For this reason a "hole" is placed in the 
flare kernel having a radius Rmin- Rmin is approximately given by X/NA. All flare 
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contributions from within Rmin are treated with partial coherent bilinear convolution, 
which is calculated within the ROL 

Division of Incoherent Flare Calculation into Long and Intermediate-Ranges 
For long-range effects, such as the case when the power law kernel is infinite in extent (y < 
2), the sector convolution must be applied on the full chip. However, it is impractical and 
not necessary to do sector convolution on all polygons on the mask. For efficiency 
reasons, the fiill mask pattern M is divided into two domains. Note that, unlike the 
subdivision of the point spread function between aberration effect and flare, the flare 
kernel is not divided into long-range and intermediate-range; instead the mask is 
partitioned. First, as depicted in Fig. 1 1 A, the long-range convolution on MI is done via a 
Fast Fourier Transform on a coarse grid 110, with each grid point 112 representing an 
average pattern density of patterns inside that particular grid point. This step must be done 
once per mask, since the pattern density in each grid point is assumed not to change 
significantly even after an optical proximity correction. Second, the near range 
convolution is done via polygon convolution inside domain M2, as shown in Fig. 1 IB. Of 
course, the domain M3, which is of the same size of M2, needs to be subtracted in order 
not to double count the contribution. Linearity and additivity are exploited in the 
convolution operator as shown below. Here the mask intensity function and the kernel are 
defined as: 

M(F) = |m(F)|'; (31) 
Kemel(F) = |/i(r)|';and, 

Flare = Kernel ® M 

= Kernel® (Ml -M3+M2) 

= Kernel ® Ml - Kernel ® M3 + Kernel ® M2 
where, ® represents a convolution operation. 

M2 is the domain defined by a predetermined distance +/- R,T,id in the x and y directions. 
Ml defines the rest of the mask. In order not to double count the polygon contribution in 
Ml and M2, polygons that lie in both domains need to be clipped at the domain 
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boundary. M3 represents the hole that is trimmed out by M2 and can be calculated easily. 
A polygon pinning algorithm is then applied. 

In Fig. 11, the total mask M is divided into Ml and M2. M2 is shown as a rectangle of 
size defined by Rmid. Mask Ml denotes the rest of the mask. Each grid point Xk on Ml has 
an average density pk. Mask M2 shows polygons 1 14 inside a range defined by +/- Rmid. 
Polygons that overlap both Ml and M2 need to be clipped and pinned on the boundary, 
like the rectangle 116 in the upper right comer. M3 is the mask portion that needs to be 
subtracted later to avoid double counting. Fig. 12 is a schematic depicting the spatial 
extent of the short range and mid-range calculation. Fig. 12A shows that usually a guard 
band is required outside the evaluation region for the short-range calculation, thus 
determine the periodicity parameter in the calculation to satisfy the pseudo-periodic 
requirement. Fig. 12 B shows that for the same evaluation box, we need to include a 
larger area for intermediate range calculation since the extent of the kernel is much 
larger. 

While the present invention has been particularly described, in conjunction with a 
specific preferred embodiment, it is evident that many alternatives, modifications and 
variations will be apparent to those skilled in the art in light of the foregoing description. 
It is therefore contemplated that the appended claims will embrace any such ahematives, 
modifications and variations as falling within the true scope and spirit of the present 
invention. 

Thus, having described the invention, what is claimed is: 
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